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Using the large hadroproduced charm sample collected in experiment E791 at Fermilab, we report 
the first directly measured constraint on the decay-width difference AF for the mass eigenstates of 
j_j ■ the D°-D° system. We obtain our result from lifetime measurements of the decays D° — > K~ir + 

and D° — > K~ K + , under the assumption of CP invariance, which implies that the CP eigenstates 
and the mass eigenstates are the same. The lifetime of D° — » K~ K + (the CP-even final state) 
is tkk = 0.410 ± 0.011 ± 0.006 ps, and the lifetime of D° -» K~ty + (an equal mixture of CP- 
odd and CP-even final states) is tkw = 0.413 ± 0.003 ± 0.004 ps. The decay-width difference is 
Ar = 2(Fkk — Fkw) = 0.04 ± 0.14 ± 0.05 ps -1 . We relate these measurements to measurements of 
mixing in the neutral D-meson system. 

PACS numbers: 14.40 Lb, 12.38.Qk 

The two mass eigenstates of the D°-D° system may have different lifetimes. In the Standard Model the domi- 
nant contribution to such a difference arises from box-diagram amplitudes involving the exchange of two W bosons. 
Suppression factors due to the GIM mechanism and heavy quark symmetries in QCD lead to the conclusion that 
these lifetime and mass differences are too small to be measured Q. Additional Standard Model contributions from 
penguin diagrams Q and long-distance effects || are also expected to be negligible. Thus, a measurable difference in 
the lifetimes [Q or masses would constitute evidence for new physics. In this letter, we report the first direct search 
for a lifetime difference in the neutral D-meson system. 

Because lifetime and mass differences lead to particle-antiparticle mixing, previous experiments have looked for 
mixing as evidence of such differences. The rate of mixing is usually characterized by the parameter r m i x , which is 
defined as the ratio of mixed to unmixed decays 
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This ratio can have contributions from both the mass difference Am and the decay-rate difference Ar between the 
two neutral D mass eigenstates according to 



_ (Am)' (AIT m 

'mix — ^ gp2 ' 

where T is the average of the decay rates of the two neutral D mass eigenstates. 

To date, particle- antiparticle mixing has only been observed in the strange- and bottom-quark sectors. Previous 
searches for D°-D° mixing have set a limit r m j X < 0.50% at 90% confidence level (CL) ||. Inserting this value into 
Eq. ||and assuming Am — 0, one obtains the limit |Ar| < 0.48 ps _1 at 90% CL ||. The measurement of Ar reported 
here is more sensitive than these previous measurements. 

If CP conservation holds in the P>°-meson system, the even and odd CP eigenstates and D® will be the mass 
eigenstates; i.e., they will have definite masses mi and m.2, and decay widths Ti and r 2 . Since the Cabibbo-suppressed 
final state K~K + is CP-even, it results from the decayf] of D\ and will have an exponential decay-rate distribution 

T KK (t) = A KK e- Tlt (3) 

where here and in Eqs. |I] and |?] below, A represents all time-independent factors. In contrast, the Cabibbo-favored 
final state K~ir + receives roughly equal contributions from D® and D\. Under the assumption of CP conservation 
and ignoring doubly-Cabibbo-suppressed terms, the decay-rate distribution of a combined sample of D° — > K~tt + 
and D° — * K + ir~ (from mixing), and separately a combined sample of D° — > K + ir~ and D° — > K~ir + , can be 
described by [@,§0 

r Xir (t) = A K ^e~ Tt cosh > (4) 

where 

r = (r! + r 2 )/2 and Ar = r!-r 2 , (5) 

and so 

Ar 11 

— = Tr - T = T KK -T K „= . (6) 

^ tkk t K 7t 

Given the current limits on Ar, Eq. |] is well- approximated by a pure exponential distribution 

r«r(t) - A K «e- vt (7) 

over the range of sensitivity of our experiment. 

If the mass states are not CP eigenstates, i.e., CP violation is present, then interpreting the measured lifetime 
difference in terms of Ti and T 2 is less straightforward. However, even in the absence of strict CP conservation, pure 
exponential distributions are still good approximations for both the KK and Kn distributions in our data sample. In 
this case, a measured difference at our level of sensitivity would constitute evidence for physics beyond the Standard 
Model, although it might not be easily related to the lifetime difference of the mass eigenstates. For convenience of 
interpretation, we hereafter assume CP invariance in the Z?°-meson system, though the validity of that assumption 
has previously been questioned in the literature || . 

The results reported here are based on data accumulated by experiment E791 in a 500-GeV/c tt~ beam during the 
1991/92 Fermilab fixed-target run. E791 was the fourth in a series of charm experiments performed in the Fermilab 
Tagged Photon Laboratory. The E791 spectrometer jl0| was an open-geometry detector with 23 planes of silicon 
microstrip detectors (6 upstream and 17 downstream of the target), 35 drift chamber planes, 10 proportional wire 
chambers (8 upstream and 2 downstream of the target), two magnets for momentum analysis, two large multicell 



1 We imply charge-conjugate decay modes throughout this letter. 
2 Reference [HI omits the factor of 2 in the cosh term. 
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threshold Cerenkov counters for charged particle identification, electromagnetic and hadronic calorimeters for elec- 
tron/hadron separation as well as for online triggering, and a fast data acquisition system that collected data at a 
rate of up to 30 Mbyte/s with a 50 /zs/event deadtime. The target consisted of a 0.52-mm platinum foil followed by 
four 1.6-mm diamond foils. Each target center was separated from the next by about 1.5 cm, allowing observation 
of charm-particle decays in air with minimal background from secondary interactions in the targets. The very loose 
transverse-energy trigger was based on the energy deposited in the calorimeters and was highly efficient for charm 
events. Over 2 x 10 10 events were recorded during a six-month period. 

The measurement of Ar is based on the observation of D° — + K~K + and K~ir + decay modes in the same detector 
with similar detector and topological biases and uncertainties, leading to a precise comparison of the lifetimes for the 
two modes. This is a companion analysis to that of Ref. |TlJ] , in which we set limits on CP violation and measured 
the branching ratios for D° — > K~K + and D° — » 7r~7r + relative to that for D a — » K~ir + . 

Selection requirements for inclusion in the D° — » K~ K + and D° — » K~ir + samples included a minimum significance 
of separation of the candidate decay vertex from the production vertex in the beam direction (> 8<r where a is the 
uncertainty in the separation of the two vertices). Candidate decay tracks were required to be inconsistent with 
coming from the D° production vertex. The net momentum of the D° candidate transverse to the line connecting the 
production and decay vertices was required to be less than 0.35 GeV/c, and the sum of the pi of the decay tracks was 
required to be greater than 0.25 (GeV/c) 2 , with p t measured relative to the direction of the candidate D°. Finally, 
to remove secondary interactions, the decay vertex had to be located at least 4er outside the target foils, where a is 
the uncertainty in the position of the decay vertex in the beam direction. 

We used particle identification from the Cerenkov counters to improve the statistical significance of the D° — * K~K + 
sample. To minimize systematic effects in calculating Ar, similar particle identification criteria were applied to the 
K~ir + sample as well. We obtained the necessary particle-identification efficiencies from a study of the D° — > K~tt + 
sample. We corrected for the Cerenkov counter efficiency by a weighting procedure whereby each two-track candidate 
decay was weighted by the inverse of the product of the two particle-identification efficiencies calculated for their 
individual p, p t , and particle type. Figure |] shows the weighted invariant mass plots for our final data sample. 
The numbers of unweighted signal events in the K~K + and K~tt + samples are approximately 3200 and 35,400, 
respectively. 

To account for the effects of our selection criteria, we define a reduced decay length for each D° candidate as the 
distance traveled by the candidate beyond that required to survive our selection criteria, and a reduced proper time 
corresponding to that reduced decay length. To the extent that the K~K + final state is a mass eigenstate of the 
neutral Z?-meson system, its proper-decay-time distribution is purely exponential. In this case, its reduced-proper- 
decay-time distribution is also purely exponential with the same characteristic decay time. In contrast, the K it final 
state is a mixture of mass eigenstates with its proper-decay-time distribution given by Eq. ^. However, since the 
difference in the decay rates of the eigenstates is small (as known from limits on r m ; x ), the difference in the amount 
of D\ and D° contributing to the D° — > K~n + decay after our proper-decay-length cut is negligible. Also, effects 
of fitting to a single exponential of reduced proper time are smaller by an order of magnitude than the systematic 
errors described below. We use reduced proper lifetime because it minimizes acceptance corrections and associated 
systematic errors. We measure the lifetimes over the same range of reduced proper time for both decay modes. 

For each bin of reduced proper time, candidates are taken from the total sample shown in Figure [l], where the signal 
region is defined as within ±2.5 Gaussian standard deviations of the Gaussian mean, and the background is taken from 
sidebands, but at the level of the cross-hatched area. We use the same D° mass and Gaussian width in each of the 
16 bins of reduced proper time, as determined from Figure separately for D° -» K~K+ (1866.7 and 11.6 MeV/c 2 ) 
and D° -» R-7T+ (1867.0 and 13.5 MeV/c 2 ). 

We correct for small acceptance differences between the K~K + and K~tt + final states by using a Pythia-based 
Monte Carlo model |Q that incorporates a full detector simulation. The simulated events are subjected to exactly 
the same selection criteria and binning as are the data. Figure ^ shows the efficiencies from the Monte Carlo study 
as functions of the reduced proper decay time for K K — + and K~tt + separately. Figure || presents the exponential 
fits to the measured reduced-proper-decay-time distributions after particle-identification weighting and acceptance 
corrections. We measure decay widths Tkk — 2.441 ± 0.068 ps^ 1 and Tkit — 2.420 ± 0.019 ps -1 ; or lifetimes 
tkk — 0.410 ± 0.011 ps and tk-k = 0.413 ± 0.003 ps, where the quoted errors are statistical only. Using Eq. [| we 
calculate Ar = 0.04 ± 0.14 ps -1 . We have also performed a joint fit to the distributions using Eqs. || and || which 
yields the same result. 

Systematic errors are studied separately for the two lifetimes and for the lifetime difference. These studies include 
the effects of varying the event selection criteria, the D° production model used in the Monte Carlo, the particle- 
identification weighting procedure, the effect of fixing the K~ K + width, and fixing the range of reduced proper 
lifetime over which the fits are done. Estimates of the systematic errors arising from these sources are summarized 
in Table |. In obtaining systematic errors on Ar from those on the lifetimes, correlation effects are included. For 
example, different production models lead to similar changes in acceptance for the K~ K + and K~tt + modes, resulting 
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in a cancellation in the systematic error in AI\ On the other hand, the uncertainties in mass resolution for the K~K + 
and K~n + samples come from fits to independent data samples and are uncorrelated, and thus their contribution 
to the uncertainty in Ar is calculated as the quadratic sum of contributions coming from the uncertainties in the 
two lifetimes. Systematic uncertainties due to the different absorption cross-sections for K and 7r are negligible, both 
because the momentum distributions are quite similar and because the targets are very thin. 

In summary, we measure tkk — 0.410 ± 0.011 ± 0.006 ps and tk-r — 0.413 ± 0.003 ± 0.004 ps. These results are 
compared with previous measurements and the world average in Table ||. We find Ar = 0.04 ± 0.14 ± 0.05 ps -1 , 
leading to a limit of —0.20 < Ar < 0.28 ps -1 at 90% CL. The value of Ar is consistent with zero and thus, at our 
level of sensitivity, is consistent with the Standard Model. This measurement gives us, independent of mixing studies, 
a limit on Ar. This 90% CL limit on the value of Ar does not saturate the 90% CL limit on r m ; x . If r m i X were at its 
limit, this would have to be due to a contribution from (Am) 2 . 

We gratefully acknowledge the assistance of the staffs of Fermilab and of all the participating institutions. This 
research was supported by the Brazilian Conselho Nacional de Desenvolvimento Cienti'fico e Tecnologico, CONACyT 
(Mexico), the Israeli Academy of Sciences and Humanities, the U.S. Department of Energy, the U.S. -Israel Binational 
Science Foundation, and the U.S. National Science Foundation. Formilab is operated by the Universities Research 
Association, Inc., under contract with the United States Department of Energy. 



[1] E. Golowich, Proceedings of the Conference on B Physics and CP Violation, Honolulu, HI (March, 1997), and references 
therein. 

[2] A. A. Petrov, Phys. Rev. D 56, 1685 (1997). 

[3] E. Golowich and A. A. Petrov, Phys. Lett. B427, 172 (1998); J. F. Donoghue, E. Golowich, B. R. Holstein, and J. 
Trampetic, Phys. Rev. D 33, 179 (1986); H. Georgi, Phys. Lett. B297, 353 (1992); T. Ohl et al, Nucl. Phys. B403, 605 
(1993). 

[4] T. Liu, Ha rvard University Report HUTP-94/E021, |hep-ph/9408330| (July, 1994); T. Liu, Princeton University Report 

HEP/95-6, |hep-ph/95084iq (August, 1995). 
[5] See for example, previous results from the E791 experiment: E. M. Aitala et al, Phys. Rev. D 57, 13 (1998); E. M. Aitala 

et al, Phys. Rev. Lett. 77, 2384 (1996). 
[6] Particle Data Group, Eur. Phys. J. C 3, 497 (1998). 
[7] I. Dunietz, Phys. Rev. D 52, 3048 (1995). 
[8] W. F. Palmer and Y. L. Wu, Phys. Lett. B350, 245 (1995). 

[9] G. Blaylock, A. Seiden and Y. Nir, Phys. Lett. B355, 555 (1995); L. Wolfenstein, Phys. Rev. Lett. 75, 2460 (1995). 
[10] J. A. Appel, Ann. Rev. Nucl. Part. Sci. 42 (1992) 367; D. J. Summers et al, Proceedings of the XXVII th Rencontre de 
Moriond, Electroweak Interactions and Unified Theories, Les Arcs, France (15-22 March, 1992) 417; E. M. Aitala et al, 
submitted to Phys. Rev. D, Fermilab-Pub-98-297-E, hep-ex/980929 (September, 1998). 
[11] E. M. Aitala et al, Phys. Lett. B421, 405 (1998). 

[12] H.-U. Bengtsson and T. Sjostrand, Comp. Phys. Comm. 46, 43 (1987). 
[13] P. Frabetti et al, Phys. Lett. B323, 459 (1994). 
[14] J. R. Raab et al, Phys. Rev. D 37, 2391 (1988). 



4 



Ktt wt'd signal = 60581 ± 353 
<Weight> = 1 .71 



2.15 2.2 




2.15 2.2 



K"K + Mass (GeV/c 2 ) 



FIG. 1. Invariant mass distribution for candidates weighted by particle identification efficiency (see text), (a) D° — > K~tt + 
decays. The solid line corresponds to a fit to a Gaussian signal distribution plus a third-order polynomial for combinatoric 
background, (b) D° — > K~ K + decays. The solid line corresponds to a fit to a Gaussian signal distribution plus an asymmetric 
Breit-Wigner (D° — > K~tt + reflection), a symmetric Breit-Wigner (D° — > K~tv + tv° reflection), and a linear combinatoric 
background. The cross-hatched areas are the estimated backgrounds beneath our signal region of width ±2.5 Gaussian standard 
deviations about the Gaussian mean. 
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FIG. 2. The K ir + (top) and K K + (bottom) efficiencies as functions of the reduced proper decay time. Note that the 
vertical scales have suppressed zeros. 
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FIG. 3. Distributions and the exponential fits for the number of D° — > K 7r + (top) and K K + (bottom) decays as a 
function of reduced proper decay time, after particle identification weighting and acceptance corrections. 
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TABLE I. Contributions to the systematic errors in lifetimes and lifetime difference. 



Systematic error in 


tk-k (ps) 


tkk (ps) 


Ar (ps- 1 ) 


Fit Range 


0.002 


0.003 


0.024 


Selection Criteria 


0.001 


0.002 


0.020 


Particle ID Weighting 


0.001 


0.003 


0.024 


MC Production Model 


0.003 


0.003 


0.000 


Fixed Width 


0.001 


0.002 


0.030 


Total 


0.004 


0.006 


0.050 



TABLE II. Comparison of our measured lifetimes with previous measurements. 





TKix (PS) 


tkk (ps) 


E791 
E687 [ 
E691 [ 
PDG [ 


Li 
5f 


0.413 ±0.003 ±0.004 
0.413 ± 0.004 ± 0.003 
0.422 ± 0.008 ± 0.010 
0.415 ±0.004 


0.410 ±0.011 ±0.006 
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